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AXIS  SYSTEM  AND  NOTATION 


The  axis  system  used  in  the  analysis  is  a  booty  axis  system,  the 
origin  of  which  is  at  the  aircraft  center  of  gravity.  The  X  axis  lies 
in  the  plane  of  symmetry,  perpendicular  to  the  rotor  shafts,  positive 
direction  forward.  The  Y  axis  is  perpendicular  to  the  plane  containing 
X  axis  and  the  rotor  shafts,  positive  to  the  right  of  the  helicopter. 
The  Z  axis  is  parallel  to  the  rotor  shafts,  positive  direction  down. 

The  symbols  and  physical  quantities  used  in  the  analysis  are  as 
follows : 


GENERAL 
(   )x 


(   ) 


D(   ) 


Bar  denotes  _=_  (   ) . 

Subscript,  following  standard  practice,  denotes  partial 

differentiation  with  respect  to  subscript  quantity; 

Dot  denotes  differentiation  with  respect  to  time; 

(')   -     d(  .  } 
(  >   "  XT  * 

Denotes  differentiation  with  respect  to  non-dimensional 
tlwj      D(      )   =  -j^.     . 


FORCES  AND  MOMENTS 


IT 
Y 

L,  N 


H 


My 


Sum  of  the  thrusts  of  the  two  rotors  (lbs.)« 

Side  force,  positive  to  right  (lbs.). 

Rolling  and  yawing  moments  of  helicopter  about  its  C.  G.; 

roll  is  positive  for  right  side  down,  and  yaw  for  nose 

right  (ft.-lbs.). 

Stabilizer  hinge  moment,  positive  for  leading  edge  right 

(ft.-lbs.). 

Rolling  moment  equation  forcing  function,  i.e.,  rolling 
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moment  due  to  lateral  cyclic  stick  displacement  (ft. -lbs.)' 
YF  Side  force  equation  forcing  function,  i.e.,  side  force  due 

to  lateral  cyclic  stick  displacement  (iba.), 
>^Cm         Coefficient  of  combined  thrust  of  front  and  rear  rotors. 

iCT  "  ,  t  R^.  l  H)2    *  r  rr  &(n  R)2   " 
Cy         Side  force  coefficient.  Cy  =  - — ngf     n^2    ' 

C,         Rolling  moment  coefficient.  C  ■  _i __  . 

1  /oTtR2^/  R)2Ji 

Q  Yawing  moment  coefficient.  Cn  =     .r . 

n  /-rr  R2(^  r)2, 


Ct         Stabilizer  lift  coefficient. 

•L'S 


M- 


Side  force  forcing  function  coefficient.  Cv  =    A*.        , , ,   • 
F  ^  XF   ^  rr  R*(.iR)id 

Cv,         Rolling  moment  forcing  function  coefficient. 
C       =     % 

Cv.         Stabilizer  hinge  moment  coefficient.  Ch  =  § _  m^    • 

a  n  .      R   r>  / 


*t  ss 


c, 


PHYSICAL  CHARACTERISTICS 


b  Number  of  blades  per  rotor,  b  =  3* 

c  Average  blade  chord,  c  =  1.5  ft. 

c8         Chord  of  vertical  stabilizer.  cs  =  k.b'J   ft, 

d^         Distance  from  aerodynamic  center  of  vertical  stabilizer  to 

hinge  line,  positive  for  a.  c.  aft  of  hinge  line,  in  chord 

lengths,  d-,  =  xac  -  XHL  =  0.25  -  *£&  . 

cs  us 

d2         Distance  from  C.  G.  of  vertical  stabilizer  to  hinge  line, 
positive  for  C.  G.  aft  of  hinge  line,  in  chord  lengths. 

<*  =  CG  ~  XHL  =  o.l»25  -  l^gL  . 
cs  cs 

h^         Distance  from  X  axis  to  spanwise  center  of  vertical  stabilizer. 
ht  =  5.0  ft. 
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Is         Mass  moment  of  inertia  of  vertical  stabilizer  about 

p 
vertical  line  through  its  hinge  point.  I8  *  3-73  slug  ft. 

Ix         Mass  moment  of  inertia  of  helicopter  about  its  X  body  axis. 

Ix  -  5^0  slug  ft.2 
Iz         Mass  moment  of  inertia  of  helicopter  about  its  Z  body  axis. 

Iz  =  80,000  slug  ft.2 

IX2        Mass  product  of  inertia  of  helicopter  for  the  X  and  Z  body 

2 


axes.  Ixz  ■  1600  slug  ft.' 
K  Tab  gear  ratio,  — -£  . 

o  6 

ks         Radius  of  gyration  of  vertical  stabilizer  about  vertical 

line  through  its  hinge  point.  k32  =  Is.  =  1.77  ft.2 

mg 

kx         Radius  of  gyration  of  helicopter  about  its  X  body  axis. 

kx2  =  hL  =  13.4  ft.2 

kz         Radius  of  gyration  of  helicopter  about  its  Z  body  axis. 

k22  =  h.  m   198.7  ft.2 
m 

X  Distance  between  rotor  shafts.  J.   =  J+2.3  ft. 

i0         Average  distance  from  C.  G.  of  helicopter  to  rotor  shafts. 

i0  »  iL  ,  21.15  ft. 

i^         Distance  from  C.  G.  of  helicopter  to  hinge  point  of  vertical 

stabilizer,  measured  parallel  to  X  axis.  X^  =  19*7  ft. 
m         Mass  of  helicopter,  m  =  J+03'7  slugs. 
mg         Mass  of  vertical  stabilizer.  mB  =2.11  slugs. 

q.  Dynamic  pressure,  q  «/|  V^  . 

2 
q^         Dynamic  pressure  at  tail  surfaces .  c^  »  /-V^. 

R  Blade  radius.     R  =22.0  ft. 

Sa  Vertical  stabilizer  area.  SQ  =  31.5  ft.2 

W  Helicopter  gross  weight.  W  ■  13,000  lbs. 

Wg  Weight  of  vertical  stabilizer.  Ws  =  67.9  lbs. 
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VELOCITIES  AND  ANGLES 

p,r        Non-diiaensionalized  roll  and  yaw  rates,  respectively. 

V  Relative  wind  velocity,  equal  approximately  to  forward 

flight  velocity  ( ft. /second) . 
oc  Angle  of  attack,  defined  as  angle  between  relative  wind 

and  a  normal  to  the  rotor  shaft  ( radians ) . 
/3  Sideslip  angle  of  helicopter,  positive  for  slip  to  right 

( radians ) . 
°6         Stabilizer  deflection  angle,  measured  between  stabilizer 

plane  of  symmetry  and  XZ  plane,  positive  for  leading  edge 

right  ( radians ) . 
$ .  f  Roll,  positive  for  right  side  down,  and  yaw,  positive  for 

nose  to  right,  angles  of  helicopter,  respectively  (radians) 
j\  Angular  velocity  of  rotor,  assumed  constant. 

fl  =  27.0  radians /second. 
AERODYNAMIC  AND  NON-DIMENSIONAL  PARAMETERS 

a  Slope  of  rotor  blade  two-dimensional  lift  curve. 

a  ■  5 '73  Per  radian. 
a8         Slope  of  vertical  stabilizer  lift  curve. 

as  =1.83  per  radian. 
h8         Non-dimensional  vertical  stabilizer  inertia  parameter. 

hc  =  ^f  C±l\    =  0.00086. 


hx         Non-dimensional  roll  inertia  parameter. 
hz         Non-dimensional  yaw  inertia  parameter. 


hx  m  I  \-%&-\     .  0.0153. 


hz  =  __i_ /_jSi).  0.227. 
hX2        Non-dimensional  product  of  inertia  parameter. 

hx2  =  ~  -£fr  -  0.«*5U. 


9. 

y[t  Tail  efficiency.   r/t  =  -1  =  0-90. 

m  Advance  ratio.  /*   =  — -^  • 

u  Helicopter  density  parameter.  /<*-   mS  .   =  6.03- 


/**, 


Air  density  (slugs/ft.3). 


cr  Rotor  solidity.   ^  =  -^  =  O.O65. 

r  Non-dimensional  time  parameter.  T  -  n  %2  R  =  0<1^  seconds 

Tj,         Tab  effectiveness.   rz   =  g-^-8 


Vertical  stabilizer  density  parameter. 


/  *o     -  So  Cr 


rB   /o  ae  cs 


10. 
SUMMARY 

An  analytical  investigation  of  the  effect  of  free  floating  vertical 
stabilizers  on  the  lateral  stability  of  the  H-21  tandem  helicopter  was 
made  for  flight  at  sea  level  and  80  knots.  The  lateral-directional  equa- 
tions of  motion  for  a  tandem  helicopter  without  free  floating  stabilizers 
were  modified  to  account  for  the  addition  of  floating  fins. 

The  basic  helicopter  stability  derivatives  were  obtained  from  the 
results  of  flight  tests  of  the  unmodified  (fixed  fin)  helicopter  by  the 
manufacturer.  The  stability  derivatives  associated  with  the  floating 
fins  were  estimated  through  a  theoretical  approach. 

The  fins  were  considered  to  be  so  mounted  on  the  aircraft  as  to  be 
mass  and  aerodynamically  unbalanced.  The  aerodynamic  unbalance  was  com- 
pensated by  a  geared  tab.  The  effect  on  lateral  stability  of  various 
combinations  of  tab  gear  ratio  and  mass  and  aerodynamic  unbalance  was 
investigated.  The  results  were  obtained  by  putting  the  modified  equations 
on  an  analog  computer.  These  results  showed  that,  although  the  floating 
freedom  of  the  vertical  fins  did  not  make  the  helicopter  stable,  damping 
of  the  unstable  rolling  mode  was  increased.  The  best  combination  of 
balance  and  gear  ratio  lengthened  the  time  to  double  amplitude  by  57$  • 
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INTRODUCTION 

GENERAL 

The  H-21  is  a  tandem-rotor  helicopter  with  three-b laded,  articulated 
rotor  heads  and  a  rotor  radius  of  twenty-two  feet.  It  flies  at  a  gross 
weight  of  approximately  thirteen  thousand  pounds,  and  is  currently  used 
by  the  armed  forces  of  several  nations  of  the  NATO  alliance.  A  drawing  of 
the  aircraft  is  shown  in  Figure  1. 

The  H-21  is  unstable,  but  pilots  feel  that,  compared  with  other  heli- 
copters, its  handling  characteristics  are  good.  Its  standard,  fixed 
vertical  tail  fins  make  it  an  excellent  test  vehicle  for  an  investigation 
of  the  possibilities  of  making  it  less  unstable,  or  even  stabilizing  it, 
through  modifications  to  these  fins.  At  the  time  of  this  writing,  the 
manufacturer  has  an  H-21  equipped  with  the  floating  fins  described  in  this 
paper,  and  is  conducting  flight  tests  with  them.  Difficulties  have  been 
encountered  in  this  flight  test  program  which  seem  to  be  due  to  aeroelas- 
ticity.  Though  this  problem  is  outside  the  scope  of  this  paper,  it  is 
mentioned  in  order  to  bring  out  the  point  that  increasing  the  size  or 
effectiveness  of  the  fins,  a  logical  extension  of  the  analysis  conducted 
here,  may  involve  practical  difficulties  which  are  not  obvious. 
DESCRIPTION  OF  FLOATING  VERTICAL  STABILIZERS 

The  proposed  free  floating  vertical  stabilizers  are  NACA  0012 
symmetrical  airfoils  of  rectangular  planform  with  a  span  of  6.75  feet  and 
a  chord  of  ^.67  feet.  They  are  mounted  on  short  outriggers  at  the  rear  of 
the  fuselage,  five  feet  above  the  X  axis,  one  on  each  side  of  the  helicopter. 
The  suspension  is  a  single  hinge  at  the  center  of  the  span,  the  axis  of  the 
hinge  being  parallel  to  the  span-wise  direction. 

The  aerodynamic  center  of  each  fin  is  approximately  at  the  25$  chord 
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point.  The  hinge  point  is  varied  between  0.28  cg  and  O.36  cg  in  this 
analysis.  The  center  of  gravity  has  been  calculated  as  0.^25  cs  . 

The  aerodynamic  unbalance  cf  the  stabilizer  is  compensated  with  a 
trailing  edge,  full  span  tab,  the  leading  edge  of  which  is  approximately 
at  its  hinge  point  which  is  at  80$  of  the  stabilizer  chord.  The  tab  is 
so  geared  that  a  deflection  of  the  stabilizer  causes  a  tab  deflection  in 
the  same  direction,  causing  a  restoring  stabilizer  hinge  moment.  The 
gear  ratio,  — -  ,  is  varied  between  1.5  and  2.5  in  this  analysis.  A 
line  drawing  of  the  fin  may  be  found  on  Page  29. 

DISCUSSION 
EQUATIONS  OF  MOTION 

The  lateral-directional  equations  of  motion  of  the  helicopter  are 
obtained  by  modifying  equations  given  in  Reference  1.  Equations  15  in 
Reference  1,  modified  to  suit  the  control-fixed  case  are: 

**■    a  <T     '  \  /U*    ex  cr  /  >-      \atr        /i*  /   -> 

Equations  1  are  now  modified  to  include  the  effects  of  the  floating 
vertical  stabilizers.  This  involves  the  inclusion  of  the  additional 
variable,  &6  ,  and  its  coefficients  in  each  equation,  and  the  formulation 
of  an  additional  equation;  that  for  stabilizer  hinge  moments. 

Equation  la  takes  on  a  term  expressing  the  side  force  due  to  stabilizer 
deflection  in  coefficient  form, 

2  Cysa6s  , 
where  Cy&g  is  developed  as  follows: 


Due  to  stabilizer  deflection,    §8   , 


Y       a  <r 


y       a^r      rrR Mt    '    5         a.  er     rr  R      Lt     5 

Taking  the  partial  derivative  of  Cy  with  respect  to  S8  ,  the  /3g  portion 
drops  out,  and 

a)     Cy  =  ^^i  -%  0  +  ii  K)jjt. 

'         r  s     acrrrRv 


This  development  has  been  for  a  single  stabilizer.  If  the  effects  of  airflow 
around  the  fuselage  are  neglected,  both  stabilizers  may  be  assumed  to  act 
identically,  and  the  total  effect  of  stabilizer  deflection  on  side  force  is 
equal  to  2  Cy  s8  ^s  * 

The  effect  of  pressure  variation  over  the  stabilizer  due  to  a  deflection 
rate,  and  that  of  the  inertia  of  the  stabilizer  and  its  angular  acceleration 
are  assumed  to  be  of  second  order  and  are  neglected. 

Equation  lb  takes  on  a  term  expressing  the  rolling  moment  due  to 
stabilizer  deflection  in  coefficient  form, 

2  C*Sb  Ss  , 
where,  by  analogy  to  the  Cy    derivation, 

3)      C   =  —  r 

Similarly,  Equation  lc  takes  on  a  term  expressing  the  yawing  moment 
due  to  stabilizer  deflection  in  coefficient  form, 

2C„SsSs  ; 

where,  by  analogy  to  the  derivations  of  Equations  2  and  3, 
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4)      0,,  =  --#  C 


Y  c 


The  effect  of  stabilizer  deflection  rate  upon  yawing  moment  may  be 
significant,  so  it  is  included  in  the  yawing  moment  equation  as 

A  general  method  of  evaluating  Cn     for  a  fixed  wing  aircraft  of  "average" 

characteristics  is  given  in  the  Appendix  of  Reference  k,   where  D'  is  a 

differential  operator,  defined  as   . ~;  .  .  This  method  is  modified  to  fit 

d(Vt) 

the  parameters  used  in  this  paper ,   and  used  as  an  approximate  method  of 

obtaining  Cn      for  the  H-21 .  A  stabilizer  deflection  rate  causes  a  change 
"D  8  8 

in  stabilizer  lift  coefficient,  A CL  .  The  resulting  yawing  moment  is 

s 

r  A  Cl*  <lt     5s  ^t 

From  Reference  k, 


o^^&'^m 


D'b, 


^m 


Substituting,  taking  the  derivative  with  respect  to  D  5 8  ,  and  multiply - 
2 


ingby— 


r            -    -A'   ^        **  ^  ^       I  ajL_        d  ^t 
C,     _  _. >   — ■  :-■*,'—■  r  a, 


'"OS*     a  a-      lt       £   m    1_<>D'6S     5  jD'i 


From  Figure  15  of  Reference  4,  and  considering  the  hinge  point  of  the 
stabilizer  to  be  at  its  leading  edge  in  using  this  figure, 

-^  =  3.15 

*  D'SS 

ifh-  =  1.50. 
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and 


There  is  a  similar  contribution  to  yawing  moment  due  to  the  tab 
deflection  rate,  K  D  6g  ; 

The  total  value  of  CnD       is,  then 
5)  C"os$~  -  fz  1*  ^^     [3.15  +  0.48K  +  0  50+o.l6K)aJ. 

The  stabilizer  hinge  moment  equation  is  now  derived.  The  tab,  being 
linked  to  the  stabilizer  in  a  manner  such  that  tab  deflection  is  a  linear 
function  of  the  stabilizer  deflection,  is  treated  as  a  modifier  in  the 
various  stabilizer  coefficients.  Tab  inertia  effects  are  neglected  as 
being  of  second  order. 

Stabilizer  hinge  moments  may  be  classified  as  aerodynamic  and  inertial. 

Hs  ■  Hsa  +  HBi  =  0 
Reference  2,  treating  the  rudder  of  a  fixed  wing  aircraft,  atates  that 
H8  -  F(/3g,  v  ,  $  ,      $,  5e,  b8), 
with  the  implied  conditions  that  the  rudder  be  mass  balanced  and  the  center 
of  pressure  of  the  rudder  be  close  to  the  X  axis.  The  vertical  stabilizer 
being  considered  is,  however,  mass  unbalanced,  warranting  the  inclusion  of 
a  function  of  roll  angle  in  the  above  relationship,  and  the  stabilizer  center 
of  pressure  is  far  enough  from  the  X  axis  to  justify  the  inclusion  of  a  func- 
tion of  roll  rate.  Including  these  two  variables  and  expanding  by  means  of 
partial  derivatives; 

b°        ^/34  ;  y  ^^  J<t  <>4> 

ih,,       ..         (3H.     ,    ,  eH-is 
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Initial  conditions  are  straight  and  level,  balanced  flight,  so  that  Hs 

and  the   A's  may  be  dropped. 

p 
Dividing  by  q^   Ss  cs  to  non-d 1  mens ionalize ,  and  multiplying  by  -=_ 

cL  ^ 

for  consistency; 

+  C,  55   +  Cw,  -  £.  . 
Treating  the  hinge  moment  due  to  inertia; 

C  =a  ^  -^  "_.     ™s    ks      &5       _     _  ?>■         w;>    k-v  I     q2   £ 

"   <x  a-      7_t  5,    cA  «ct-     jr   5,    c,    ra 


s 


CL-  -  -I—  hs  0Z  hs 


The  stabilizer  hinge  moment  equation  is  now 
6)    Ch/8>/i9  -  (lkw  +  C^  D)  D^  ♦  <c^  *  Chw  d  *  c^  D2-)  * 

The  equations  of  motion  for  the  helicopter  with  floating  vertical 
tail  fins  are: 

7b)     C,   ^s  +  (^r  4-  W  O;  W  f  (  &J1  d  -  >vdM*i-2  ^   S5  =  0 
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DEFINING  STABILIZER  HINGE  MOMENT  DERIVATIVES 

The  stabilizer  hinge  moment  derivatives  in  Equation  7d  must  now  be 
further  defined. 

Ch*    is  the  result  of  the  aerodynamic  unbalance  of  the  stabilizer.  A 
sideslip  angle,  /%B  9   it  equivalent  to  an  angle  of  attack  on  the  vertical 
stabilizer,  and 

H  -  as/3s  qt  Ss  (dx  oa) 

U  «  2  HJ*M  31  Sg  dl  0<^  «  2  ^  aS  /*\  A 
h  a  <r  q  s8  cg  (t     a   ^     1  /  s 

8>  *V     27?t^\ 


C^_   is  a  result  of  the  aerodynamic  unbalance  of  the  stabilizer  and 
the  effective  change  in  vertical  stabilizer  angle  of  attack  due  to  a  yawing 

velocity. 

H  =  a8  ^  qt  8.  (di  oB) 

vt  —7—  r    y 

h  =  -a«Att  sg  di  °f  pr 

Oh^^i^y  S§dT  o,  Vp^        ~       aj^dlADH, 
n  ao-rVqSaog  "  %  a  <r   yvX 

9)  0wr-2^il^ilA 

y'  ^  ft  acr   ^    T 

C^-2   is  a  result  of  the  mass  unbalance  of  the  stabilizer  and  its 
linear  acceleration  with  the  helicopter's  acceleration  in  yaw. 

H*  -  %  ^  V  (d2  •.)  =  -  2sJlJ2^  c2^ 
n  a  a-         q  sB  c8  a^JV7    ^        *  w  * 

10)  oh.=    -JL^A^d2 
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Ch   is  a  result  of  mass  unbalance  of  the  stabilizer.  At  any  roll 
angle  there  is  a  component  of  the  stabilizer  weight  acting  to  create  a 
hinge  moment. 

H  =  -  vvs*  dz   cs 

C    =  2/"2   Ws  dz  cs    x  __  _   2  mZ   Ms  3  <*z    is    . 4sVZ     n       S  <*Z   cs    a 


4SZ     n.      7  4^ 


Ch^.  is  a  result  of  aerodynamic  unbalance  of  the  stabilizer  and  the 
effective  change  in  angle  of  attack  of  the  stabilizer  due  to  a  rolling 
velocity.  By  analogy  to  the  derivation  of  C"hna>  > 

'         Id*         a.  cr   />*  Ja 

^hrfTt)   is  a  result  of  mass  unbalance  of  the  stabilizer  and  is  analogous 


to 


Ch^tj. 


Cjjg   is  a  result  of  the  aerodynamic  unbalance  of  the  vertical  stabilizer 
and  the  sideward  lift  resulting  from  a  deflection  (equivalent  to  a  change  in 
angle  of  attack) .  Here  the  effect  of  the  tab  must  be  included,  both  as  it 
affects  lift  and  as  it  affects  the  hinge  moment.  Considering  first  the  hinge 
moment  due  to  the  lift  of  the  deflected  stabilizer: 
H«  -  *a  U3  -/3S  +  ft   K5a)  ^tSsd,  cs 

The  additional  hinge  moment  due  to  tab  deflection  is 


'he 


^    a  j-  jj  &t 
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The  quantity,   J"1  is  evaluated  using  a  relationship  derived  in  Reference  3: 

~~  =  -   1.2606  K*   Kb  Kt  Kc  \f 

O   Ox. 

As  evaluated  in  the  CALCULATIONS  section  of  this  paper, 


i  st 


The  total  value  of  Cv,,.   is,  then 

6s 


Ik)  r,Stf  =  -Zifc  2±£l    d,    Cl+nK)  -  0.66%^  <. 


Cy,    is  a  result  of  the  change  in  the  pressure  pattern  over  the 
stabilizer  surface  with  a  stabilizer  deflection  rate.  From  Reference  k, 

r        -     { *  ^  )     +  i±S±)    (*c«  ) 
Figure  l6  of  Reference  U  gives 


2,33 


ws7A 

(1  C(0  -0  38 


and 


^r-^tx^35*0-36^ 


There  is  a  similar  contribution  to  Cn_   due  to  tab  deflection  rate, 


K  D  Ss  . 

r. 


Ch  ^  =  -  ^  J-^  fl  (0.  8  K*  O.  \A   Kas) 


Tlie  total  value  of  C^  is,  then 

"5  s 


15>  c^Ss  "  -  arJ^xL2-35  +  °-8  *  •■<  038  *  ° '«  K)  aJ 
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NUMERICAL  EVALUATION 

In  order  to  obtain  analog  computer  solutions  of  the  helicopter  lateral 
modes,  all  coefficients  are  numerically  evaluated.  The  stabilizer-fixed, 
controls -fixed,  helicopter  stability  derivatives  are  evaluated  using  flight 
test  data  furnished  by  the  manufacturer.  The  coefficients  of  the  stabilizer 
hinge  moment  equation  and  the  stabilizer  deflection  portions  of  the  heli- 
copter equations  are  evaluated  by  using  the  expressions  derived  in  the  preceding 
sections  of  this  paper. 

STABILIZER-FIXED,  CONTROLS -FIXED  STABILITY  DERIVATIVES 
Test  conditions: 

Altitude:  Standard  sea  level 

V  =  80  knots 

_n_  =  27.0  radians  /second 

Manufacturer's  flight  test  data: 

J  A        m   -13,000  lbs. /radian 
s 

Y^    -0 

Zl       i  W  =  13,000  lb. 

0£    =0 
o 

Y<j>  =  0 

L/5S  -  -22,500  ft.  lb. /radian 

L*;  ■  -  3>687  ft.  lb. /rad. /second 

L$,  =  -  6,600  ft.  lb. /rad. /second 

N^  -  20,000  ft.  lb. /radian 

N^  =»  -14,300  ft.  lb. /rad. /second 

N<j,  ■   2,200  ft.  lb. /rad. /second 
Non-dimensional  stability  derivatives: 

CY  ■  -0.05^6 
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\ 

=  0 

SCT 

-  0.0546 

cY 

=  0 

%. 

=-0.002245 

B'r 

=  -0.0103 

% 

«  -0.01882 

Cn  =  0.002245 

cn  =  -0.0398 

1  r 

Cn_  *  0.006l4 


STABILIZER-FREE  DERIVATIVES 

Cx         =  0.00475  +  O.OOI89  K 


o 


3 


C »    ■  0.000562  +  0.000224  K 
Cn*   a  -0.00221  -  O.OOO885  K 
CnD   ■  -0.00028  -  O.OOOO368  K 

S 

Ch^     -  0.459  dj_ 
3    PhDv     ■  "°-3575  d2 
Chj2^  -  -0.479  d2 
Ch<{>      ■  -0.0276  d2 
Chj^     ■     0.0906  d-L 
%£<>  =     0.1217  d2 

ChR       ■  (-0.458  -  O.I83  K)   di    -  O.I83  K 

*$  "*• 

Ch         -  -0.0334  -  0.0116  K 
DSS 


OTHER  NON-DIMENSIONAL  COEFFICIENTS 
=  1.22 

a  0.0102  9 


2^a 


a.<T 


A  Hs 


nt 
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NUMERICAL  EQUATIONS  FOR  FLIGHT  AT  SEA  LEVEL  AND  80  KNOTS 

16a)  (-0.0546-1.2  2  0)/^-  1.22  DY  +0.0546  4>  f  (o.OOSS  +  0.OO3  70  K)ss  =  0 

16b)         -0.002245/5,  +(-0.0O|96|  +  0.004540)0*  +-(-0003510  -  0.0/ J"  3  Dz)$ 

■f  (O.OOM24  +  O.O0O44&K)  6S  =  O 

16c)  0.  0020 /3S +6o.  007&- 0.22  7  D)D^»  +  (O.0OII7  D  +  0.00*54  Dz)^ 

+  [-0,00442-0.00l7  7K~(O.OOO56  +  0. 0000736  K)  d]  Ss  =  O 

16*0  0459  d,/Ss+ (-0.3^75-^,-0.47940)0^ +  (-0.027(5  41  +  0.0306,^,0 +  0.1217^^)4 

-f[(-0.45"8-OJ83K)d,  -O.I63  K  +  (-O.OJ3<i  -  0.O|  16  K)  D  -  0.0/  029  Dl)  S5  =  O 

FORCING  FUNCTIONS 

To  obtain  results  for  different  combinations  of  gear  ratio  and  sus- 
pension point  which  may  be  compared  in  a  straightforward  manner,   a  control 
input  is  used  to  provide  forcing  functions  in  the  computer  setup.     The 
control  input  used  in  this  analysis  is  a  one  inch  right  displacement  of 
the  cyclic  stick,   displace  and  hold. 

Manufacturer's  data  indicate  that,  at  80  knots  and  sea  level,   the 
rolling  moment  produced  by  lateral  cyclic  oontrol  is  119,880  ft.  lb.  per 
radian.     The  control  linkage  gives  0.015  radians  per  inch  of  stick  deflec- 
tion.    Therefore,  the  applied  rolling  moment  for  one  inch  of  lateral  cyclic 
stick  displacement  i6  1,798  ft.  lb.     Non- dimensional! zing  and  multiplying 


2 

by  g-p  for  consistency, 


Cm  «  0.0001798. 


This  becomes  the  forcing  function  on  the  right  side  of  Equation  l6b,  with 
a  minus  sign. 

The  manufacturer  also  indicates  that  right  cyclic  stick  displacement 
produces  a  side  force.  This  is  the  gross  weight  times  the  swash  plate 

displacement  in  radians,  or  0.01 5  W  per  inch  of  lateral  cyclic  stick  dis- 
placement. Non-dimensionalizing  and  multiplying  by  z       for  consistency, 

a  c 
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Cy       =  0.0008l8. 
XF 

This  becomes  the  forcing  function  on  the  right  side  of  Equation  l6a,  with 
a  minus  sign. 

TABULATION  OF  COEFFICIENTS 

Ail  coefficients  in  Equations  16  which  contain  d-^  and  cU  change  in 
value  during  the  investigation,  as  the  suspension  point  is  changed.  Simi- 
larly, coefficients  containing  K  change  when  the  gear  ratio  is  changed.  In 
order  to  make  the  equations  more  adaptable  to  computer  solution,  they  are 
rewritten  as 

17a)         -  1.22  D/S6  =  0.0546*  /"g  +    1.22  D^  -  0.0546**  -  Oj*  ,  -  0.000818 

17b)         -  0.0155  D2<t>  =     0.002245/3     -  0.00454  D2*  +  G.001961  D^ 

+  0.00551  D<t>  -  C2  iB  -  0.0001798 
17o)         -  0.227  D2V~  -  0.0020  /}B  +    0.0076  D*  -  0.00454  D24>  -  0.00117  D* 

+  °5  D  §a  -hC4  8a 
17d)         -  0.1985  D2  ha  =    -  G=/39  +■  05  D2^  +  Cy  DV  -  Cq  D24>  -  0o  D<*> 

+  c10*  ^cn  Dbs  +  <tl2  ss 

where 

Ci  ■  0.00950  +  0.00378  K 

C2  -  0.001124  +  0.000448  K 

C3  ■  O.OOO56  +  O.OOOO736  K 

C4  =  0.004te  +  0.00177  K 

C5  =  0.459  dx 

C6  =  0.479  d2 

CT  =  0.3575  di 

C8  »  0.1217  d2 

C9  =  0.0906  dj, 

C10  =  0.0276  d2 

Cll  ■  0.0334  +  0.0116  K 
C12   =  (0.458  +  O.I83  K)   dx  +  O.I83  K 
and  the  values  of  these  constants  are  tabulated  in  Table  I. 
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COMPUTER  SOLUTION 

Equations  17  are  now  set  up  on  the  L3  GEDA.  analog  computer.  The 
computer  diagram  and  potentiometer  settings  are  shown  in  Figure  2.  A 
photograph  of  the  computer  is  shown  in  Figure  3,  though  it  is  to  be  noted 
that  the  problem  set  up  on  the  problem  board  in  the  figure  is  not  the 
one  treated  here.  Four  computer  runs  are  made  for  each  gear  ratio,  each 
of  the  four  being  for  a  different  vertical  stabilizer  suspension  point. 
The  quantities  4> ,  D<fc  ,  D  ^,  and  6   are  recorded  on  the  Sanborn  recorder 
(left  foreground  of  Figure  3).  The  recorded  computer  results  are  shown 
in  Figures  4-8. 

RESULTS 
ROLL  AHGLE 

The  results  of  this  analysis  can  best  be  interpreted  by  observing 
the  roll  angle  traces  from  the  Sanborn  recorder.  Figure  4  shows  that  the 
rolling  motion  of  the  helicopter  with  fixed  fins  is  a  straight,  oscilla- 
tory divergence,  with  a  period  of  approximately  7*3  seconds  and  time  to 
double  amplitude  of  approximately  7  seconds .  Figures  5-8  show  that  the 
major  result  of  adding  the  floating  fin  degree  of  freedom  i6  to  increase 
damping.  A  secondary  result  is  a  slight  reduction  in  period.  A  new  long 
period  mode  has  also  been  introduced  into  the  roll  motion.  A  separate 
computer  run  was  made,  recording  only  roll  angle,  with  the  record  sealed 
down  to  show  up  this  long  period  mode  to  best  advantage.  The  computer 
overloaded  after  seven  of  the  seven-second  oscillations,  but  enough  of  a 
record  was  obtained  to  show  that  the  new  mode  is  a  pure  convergence,  reach- 
ing steady  state  in  approximately  four  cycles  of  the  seven-second  mode.  The 
seven-second  mode  continued  to  oscillate  divergently  about  the  steady-state 
value  of  the  pure  convergent  mode. 


25. 

Figures  9  and  10  and  Table  II  snow  the  variation  in  damping  with 
suspension  point  and  with  tab  ge^r  ratio.   It  is  evident  from  these  figures 
that  the  optimum  combination  of  aerodynamic  unbalance  and  gear  ratio  is  not 
in  the  range  examined.  An  additional  series  of  computer  runs  (records  of 
which  are  not  shown  here)  shows  that  even  further  improvement  of  the  damping 
is  effected  by  moving  the  suspension  point  farther  forward  (bringing  the 
absolute  value  of  dj  closer  to  zero),  and  by  decreasing  the  gear  ratio.   Un- 
fortunately, a  crowded  computer  schedule  and  other  factors  combined  to 
prevent  detailed  analysis  of  the  additional  computer  runs  or  continuation 
of  the  investigation  until  the  optimum  combination  of  suspension  point  and 
gear  ratio  were  found. 

In  the  area  investigated,  the  best  combination  of  suspension  point 
(d^=  -  0.03)  and  gear  ratio  (K=*1.5)  produced  an  increase  in  time  to  double 
amplitude  for  the  roll  angle  oscillation  of  approximately  ^7$ 

Figure  11  shows  the  variation  in  period  of  the  roll  angle  oscillation 
with  suspension  point  and  gear  ratio.   It  is  seen  that,  at  K  equal  to  1.5* 
and  dj_  equal  to  -0.03  (the  optimum  damping  combination  of  those  tested),  the 
period  is  reduced  by  approximately  2.1;$.   This  is  not  significant,  and  is 
smaller  than  the  limits  of  accuracy  of  the  Sanborn  recorder  record  and 
methods  for  measurement  of  quantities  recorded  on  it. 
ROLL  RATE 

The  effect  of  the  addition  of  the  floating  degree  of  freedom  to  the 
helicopter  on  the  roll  rate  modes  can  be  seen  in  Figures  k   -  6  to  be  mainly 
that  of  additional  damping,  plus  a  small  decrease  in  period,  as  in  the  roll 
angle  motions. 

Figure  12  shows  again  that  the  optimum  combination  of  suspension  point 
and  gear  ratio  was  not  included  in  the  test  runs . 

A  comparison  of  Tables  II  and  III  indicates  that,  within  the  limits 

of  accuracy  of  the  record  and  record  analysis,  the  effect  upon  the  period 
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of  the  roll  rate  is  the  same  as  upon  that  of  the  roll  angle. 
YAW  RATE 

The  yaw  rate  records  in  Figures  4-8  were  made  in  order  to  disclose 
unusual  variations  in  the  yaw  rate,  if  any.  The  problem  solution  did  not 
run  long  enough  to  permit  a  more  or  less  quantitative  analysis  of  this 
motion.   Qualitatively,  it  can  he  seen  that  the  addition  of  floating  fins 
gives  a  higher  mean  yaw  rate  and  better  damping  of  the  oscillations,  while 
shortening  the  period  slightly. 
STABILIZER  DEFLECTION 

As  in  the  case  of  yaw  rate,  the  stabilizer  deflection  was  recorded  in 
an  efiort  to  locate  any  irregularities  which  might  show  up  in  this  degree 
of  freedom.   The  only  irregularity  to  be  noted  is  an  expected  one;  the 
heavily  damped,  very  short  period  oscillation  which  is  seen  at  the  beginning 
of  each  run.   Although  adjustment  of  the  Sanborn  recorder  precludes  any 
positive  statements  about  the  magnitudes  of  the  stabilizer  deflection,  it 
is  safe  to  say  that,  over  the  time  interval  studied,  they  are  not  very  large, 
and  appear  to  tend  to  remain  well  within  the  range  of  mechanical  possibility. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  addition  of  a  free  floating  degree  of  freedom  to  the  vertical 
stabilizers  01  the  H-21  helicopter  theoretically  produces  more  favorable 
damping  characteristics  of  the  lateral  oscillations.   Varying  the  aerodynamic 
unbalance  of  the  stabilizers  varies  the  amount  of  improvement  in  damping  at- 
tainable, as  does  varying  the  gear  ratio  of  the  balancing  tab. 

Though  the  combinations  of  aerodynamic  unbalance  and  balancing  tab  gear 
ratio  tested  do  not  contribute  enough  to  stability  to  make  the  damping  posi- 
tive, the  use  of  the  optimum  combination  in  conjunction  with  other  stability- 
improving  devices  could  conceivably  bring  about  lateral  dynamic  stability. 
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In  the  atsence  of  other  stability-improving  devices,  it  is  doubtful 
that  any  combination  of  aerodynamic  unbalance  and  balancing  tab  gear  ratio 
could  bring  about  positive  damping  of  the  lateral  dynamic  modes.   However, 
further  investigation  of  the  areas  not  investigated  is  recommended,  in 
order  that  the  optimum  combination  might  be  found.   The  root -locus  method 
could  be  used  here  to  great  advantage. 

The  analysis  was  conducted  for  one  flight  speed,  and  its  applicability 
is  thereby  limited.   It  is  recommended  that  several  other  flight  speeds 
also  be  investigated  in  an  effort  to  determine  whether  or  not  the  floating 
vertical  stabilizer  has  a  negative  effect  on  lateral  stability  at  other 
speeds . 

The  low  aspect  ratio  and  small  area  of  the  stabilizers  considered  lead 
to  a  recommendation  that  the  same  problem  be  investigated  using  stabilizers 
of  increased  span,  thus  increasing  area  and  aspect  ratio  and  making  the 
stabilizers  more  effective. 
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CALCULATIONS 


2<?. 


x?     - 


— /c.-- »K- 


7i 


L 


<:*  -sh 


/-c;.  — 


t 


K4>=    1.3  -   0.026  *  <t>=i£>° 

K$  =    I. "5  ~  0.39  =  0.9 


Kt  = 


fe(S)-»'° 


Kc- 


07  +   0.51  £f 


fe) 


Note:   Drawing  is  approximately  to  scale.  Measurements  made  on  actual 
blueprints. 
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c,  =  1.3  inches  6  =  i?6 .0  inches 

cf  =  (1  -  0.25  +  d1)  56  =  (0.75  +  dj_)  56 

KQ  =  0.326  (0.75  +  dj_  )"°'T  +  0.51  (0.75  +  dx) 

Kb  -  1  -  O.85  [(^)2  "  (sfj)2  j       %  ■  (0.25  -  dj_)  56, 

t^  =  6.8  inches 
Evaluating  the  product  of  K^  and  Kc  with  the  hinge  point  of  the  stabilizer 
at  its  extremes, 

at  d±   =  -0.03,  KbKc=  O.687 
and  at  d±   =  -0.09,  K^K  =  0.60k. 
For  simplicity,  the  mean  will  be  used  for  all  bilge  points,  i.e., 

K^Kc   =  0.646. 
Then, 

C^  =  (-1.2606)   (0.9)   (0.61+6)  (0.90)  =  -0.66. 


a  st 

N0N-DIMENSI0N*LIZIN0  MANUFACTURER'S   DATA 
2  Y/3  s 

CY/3S   =  a7"      />^R^(-rt_R)2 

r^=  (5.,75K0.065)  =  5'36 

/ot^R2   =   (0.00238)(1520)    =  3.62 
(/ZR)2   =    [(27.0)(22.0)]2   =  353,000 

r^  r2(-/ir)2  =  1,228,000 

~C^s  ■  5-36  i^SMy    .     .  o.0^6 

-      .      2     W _      ,  13,000 

1CT  =  a7  TtT^AR?       =  5*36  (3-62)(353,000)   =  °'°546 

=  .  _2 Ms 

/>  TTR2(/iR)?£    =   (l,228,000)(l+2.3)    =   5^.05  x  106 


%s  =  5-365^oVFio^=-0-0022^ 

*  _  -OR  _2 L» 
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m  _   (27)(22)    _ 

Jo   ~     21.15        "  ^"^ 

cXr  .  (29.15X5.36)  5ll.o4  x  106  =  -0.0103 

c-tp-^rf  ir-  ^^i  =  (29.15x5.36)  5orno*  =  -°-01882 

~C»r  -  ^  £  ^rr^B)^     -  (29.15X5-36)  ^Tlo^  "  -°-^8 
S^^  ^r^Ar)^     ■  (29-15X5.36)  ^.g^  10b  -  0.00614 
EVALUATING  DERIVED,   STABILIZER-FREE  DERIVATIVES 

uc-    Sc     So 


'8 


a#  ~  (5.73)(0.065J   "  k'92 


Ss^  -  :l-5 


_^  =  --   =  0.0207 


V 
S  ■  7ITR- 


For  V  =  135  ft. /sec,       Jl  =  27. 0  rad./sec, 

/*  =   (27)(22)   =  °'2275^2   =  °'0518 
o     a.       S 

^    iT^frfe-  nt  =  (o.o5i8)(i+.92)(o.o207)(o.90)  =  0.001+75 

From  Figure  5-33  of  Reference  2, 

^2   =  0.40. 
CYc         =   (0.00l+75)(l  +  0.4  K)    =  0.00475   +  0.00189  K 

^S  =  T  ^«s     =  I+2~f  ^Ys     =  0-°°05o2   +  0.000224  K 

~C»*s        =  "  ~T  eY6fi   =   "  33  5YSg   =   -  0-00221   -  0.000885  K 

CnD5s     =-^^t  ^/m  t    [3.15   +  0.48  K  +  (1.5   +  0.16  K)   a^j 


52. 
C^  .  -  (0.2275X8.68)  (0.90)  <0-^8}jfey7)  £        J 

Cn     =  -0.00280  -  O.OOOO368  K 

^Ms   =  27?tZaT^-  dl  =  (2)(0.90)(J+.92)(0.0518)  d±   =  0.459  *i 

\fa    m  -  °-479  d2 

%      -   -  ~~  ~ §^  d2  "  "  (10*720)  I0*0*1*)  ^-°5)  {Fillz&'6J)te 
\.        =   -  0.0276  d2 

^    *    27tf>^^»- 51^-^  0.3575  ^.0.0906^ 

<W  ^  V  d2  »  "  |J  V'v  -   iM  °'479  d2  =  °'1217  d2 

X   -     -  27lt  £&  dl  ^    +  T2  K)  -  0.66  %4  K 

°h<0   a     -  (!•»)  (^92)  (0.0517)  (1  +  0.4  K)  dx  -  (3.54)  (0.0517)  K 

Ohg       *  (-  0^58  -  0.185  K)  dx  -  0.183  K 

^  *  "  ^5^V(2l!l^  L2^5  +  0.8  K  +  (0.38  v  0.14  K)  (1.83)] 

ChDs   =  -  0.01098^2.235  +  0.8  K  4-  0.695  +- 0.256  K]  =  -  0.0354  -  0.0116  K 


'3 


NON-DBKtf  SIGNALIZING  FORCING  FUNCTIONS 


Mp  -  (119,880  Ft.  Lb./kadian)  x  (Control  Displacement  in  Inches) 
x  (0.015  Radians  per  Inch  of  Displacement) 
Control  Displacement  —  1  Inch 
Mp=  1,798  Ft.  Lb. 

°m_  -  ~r    J%, — rrs  -  5.56    ,  hW       =  0.0001798 
"f     *^ /*rr  r2(ar)2j  54.05  x  106 

YF  -  (W)  (0.015)=  (13,000)  (0.015)  ~  195  Lb. 


55. 


-5e ■  « iS5 


TABLE  I 


VALUES  OF  COEFFICIENTS  CONTAINING  K,  dj_,  dg 


3-?. 


COEFFICIENT 


Cl 


C2 


C3 


Cl* 


C5 


C6 


C7 


C8 


C9 


ClO 


K 


1.5 

2.0 
2.5 

1.5 

2.0 
2.5 

1-5 
2.0 

2.5 
1.5 

2.0 
2.5 


-0.03 
-0.05 
-0.07 
-0.09 

-0.03 

-0.05 
-0.07 
-0.09 

-0.03 
-0.05 
-0.07 
-0.09 

-0.03 
-0.05 
-0.07 
-0.09 

-0.03 
-0.05 
-0.07 
-0.09 

-0.03 
-0.05 
-0.07 
-0.09 


0.U5 
0.125 

0.105 

O.O85 


0.ll*5 
O.125 
0.105 
O.085 


0.11*5 
O.125 
0.105 
O.O85 


VALUE 

0.01518 
0.01703 
O.OI89 

0.001797 

0.00202 

0.00226 

O.OOO67I 
0.000706 
0.00071+5 

0.00706 
0.00796 
O.OO885 

-O.OI38 
-0.022";? 
-0.0321 
-0.01*13 

O.O69I* 

0.0598 
0.0501* 

0.01*07 

-0.0107 
-O.OI785 
-O.O25O 
-0.0322 

O.OI767 
0.01521 
0.0128 
0.01036 

-O.OO272 
-0.001+52 
-O.OO63I* 
-O.OOS16 

0.001*00 
0.003 1+6 
0.00291 
0.00235 


35. 


TABLE  I  (CONTINUED) 
VALUES  OF  COEFFICIENTS  CONTAINING  K,  d-^  d2 


COEFFICIENT 

K 

— —  ■— — ■■              

dl 

dg 

VALUE 

Cll 

1.5 
2.0 

2.5 

0.0509 
O.O567 
0.0624 

Cl2 

1-5 

-0.03 
-0.05 
-0.07 
-0.09 

0.253 
0.238 
0.224 
0.209 

2.0 

-0.03 
-0.05 
-0.07 
-0.09 

0.3142 
0.325 
0.309 
0.292 



2.5 

-0.03 
-0.05 
-0.07 
-0.09 

0.431 
0.413 
0.395 
0.375 

3  b. 
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FtG. 2 
\-\-Z\   Geda  Computer   Diagram 

FlO*T|NJO   VCffTICAL    Tail.   Fins  56a  Level 

Controls    Fixed  60  Knots 

On£"  Inch    Lateral.    Cyclic    Displace    and    Hold   Control     Input 

FOUATlONS  : 

Y:  -  1.  a  2  D/\  =  0.0346  ft    +  1.22   Of  -  0.0346  4>  -  C,  55   -  0.  0O08I  8 

L:  -0.0153  U1*=  0.002245/3,    -000454  D^    +-O-0OI96 1  O^  *  0.  00351  D4>  -  C,  6,  -OOO0I798 

IS*  -0.227   0UV  =-O.0020  A  +0.00/6  D*    -  0.  00454  D1*    -  0.  001  I  7  D*    +  C3  D64  +C4  6, 

H-  -O.I985DlSs=-C,/3t    +C,CiV*C70*-C,Dl*-CJ0**C„t+Cl,Oo,i-C,1*i 
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